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Abstract. We present a nonextensive version of the QCD-based Nambu - Jona- 
Lasinio (NJL) model of a many-body field theory describing the behavior of strongly 
interacting matter. It is based on the nonextensive generalization of Boltzmann-Gibbs 
(BG) statistical mechanics used in the NJL model, which was taken in the form 
proposed by Tsallis characterized by a dimensionless nonextensivity parameter q (for 
q — > 1 one recovers the usual BG case). This new phenomenological parameter accounts 
summarily for all possible effects resulting in a departure from the conditions required 
by application of the BG approach, and allows for a simple phenomenological check 
of the sensitivity of the usual NJL model to such effects (in particular to fluctuations 
of temperature and correlations in a system of quarks). As an example, we discuss 
the sensitivity of such a g-NJL model to the departures from the NJL form, both for 
q > 1 and q < 1 cases, for such observables as the temperature dependencies of chiral 
symmetry restoration, masses of vr and a mesons and characteristic features of spinodal 
decomposition. 



PACS numbers: 21.65. Qr; 25.75.Gz; 05.90.-t-m; 11.30.Rd 



Nonextensive effects in the Namhu-Jona-Lasinio model of QCD 



2 



1. Introduction 

In all standard studies of high energy collisions and properties of nuclear matter one 
often uses a statistical approach based on Boltzmann-Gibbs (BG) statistics. However, 
such an approach is, strictly speaking, only correct when the corresponding heat bath 
is homogeneous and infinite. These conditions are not met in realistic situations in 
which one encounters some inherent problems arising, for example, from the smallness 
of the collision system and from its rapid evolution. These, among other things, render 
the spatial configuration of the system being far from uniform and prevent global 
equilibrium from to be established (cf. [1] and references therein). As a result, some 
quantities become non extensive and develop power-law tailed rather than exponential 
distributions. 

However, one can still use a reasonably simple statistical approach provided it will 
be based on a nonextensive extension of the usual BG statistics known as g-statistics (cf. 
[21 El m [5] and references therein) . The new phenomeno logical nonextensivity parameter 
q occurring here accounts for all possible dynamical factors violating assumptions of the 
usual BG statistics which is recovered in the limit of g ^ 1. Because it enters into the 
respective formulas of the particular dynamical model used for a given investigation, it 
allows for the phenomenological check of the stability of the model considered against 
possible deviations from the BG approacljji. 

In what concerns possible physical interpretations of the parameter q, the most 
popular one is that for the g > 1 case, g — 1 is a measure of intrinsic fluctuations 
of the temperature in the system considered [HI [151 H] whereas g < 1 is usually 
attributed to some specific correlations limiting the available phase space [16] or to the 
possible fractality of the allowed phase space [IT] (still other possible interpretations 
were considered in [I3]). As a result of application of the g-statistics, one gets a 
characteristic power-law distribution in energy-momentum [^ and specific g-versions 
of the Fermi-Dirac (FD) distribution [T8l|19] (see also [20119]). 

Notice that in g-statistics we do not specify what is the dynamical origin of these 
intrinsic fluctuations or specific correlations. It is expected that every piece of a new 
dynamical knowledge accumulated during systematic studies of the respective processes 
substantially lowers the values of the parameter |g — 1| needed to fit experimental data. 
This was confirmed in [7| when investigating transverse momenta distributions in heavy 
ion collisions, namely the gradual accounting for the intrinsic dynamical fiuctuations 
in the hadronizing system by switching from pure statistical approach to the modified 
Hagedorn formula including temperature fiuctuations [/7j , resulted in sizeable decreasing 

I Applications of the nonextensive approach are numerous and cover all branches of physics [3J. For 
those in high energy multiparticle production processes see [4j [6l [7] whereas O [3 [9] deal with different 
aspects of nuclear and quark matter. The nonextensive framework can also be derived from a special 
treatment of kinetic theory investigating complex systems in their nonequilibrium stationary states 
[SI [TOl [11] . In [12] one can find updated discussion of power-law tailed distributions emerging from 
an maximum entropy principle. Some examples of more specialized topics can be found in [13' and 
references therein. 



Nonextensive effects in the Namhu-Jona-Lasinio model of QCD 



3 



of g — 1. This means, therefore, that when one reaches in such a procedure the value 
g = 1, it should signal that all dynamical effects spoiling the initially assumed BG 
approach have already been successfully accounted for. 

Recently, g-statistics has been applied to the Walecka many-body field theory [2T] 
(known as quantum hadrodynamics or QHD-I) [22j. It resulted, among other things, in 
the enhancement of the scalar and vector meson fields in nuclear matter, in diminishing 
of the nucleon effective mass and in hardening of the nuclear equation of state (only the 
g > 1 case was considered there). 

In this paper we shall present a nonextensive version of the QCD-based Nambu 
- Jona-Lasinio (NJL) model of a many-body field theory describing the behavior of 
strongly interacting matter by accordingly modifying the NJL model recently presented 
in |23]. This means than that, unlike in [22], we shall work on the quark rather than 
the hadronic level. Also, unlike in [22], we shall consider both the g > 1 and g < 1 
cases. This will allow us to discuss the g-dependence of the chiral phase transition in 
dense quark matter, in particular the quark condensates and the effective quark masses. 
Their influence on the masses of vr and cr mesons and on the spinodal decomposition 
will be also presented. In the present work we shall limit ourselves to investigation of 
the response of these two elements to the departure from the usual BG approach and 
confront the obtained results with a possible dynamical explanations 

Our paper is organized as follows: In Section [2] we provide a short reminder of the 
basic features of the NJL model used in [23] . Section [3] contains a formulation of the 
nonextensive version of the NJL model (the g-NJL) whereas our results are presented in 
Section H] where we discuss the influence of nonextensive statistics on chiral symmetry 
restoration in g-NJL (Section 14. ip and the g version of the spinodal decomposition 
(Section 14.21) . We close with a summary and conclusions in Section [51 



2. Basic elements of the NJL model 



The SU(3) NJL model with U(1)a symmetry breaking was first formulated and discussed 
in [2^ and the first bosonized version of the NJL that obeys all strictures of chiral 
symmetry was discussed in [25]. We start with recollecting some basic formulas 
concerning the NJL model used in [23] (see also [261 1211 129]). They used the 
usual lagrangian of the NJL model, invariant (except of the current quarks mass term) 
under the chiral SUL(3)®SUij(3) transformations (described by coupling constant gs) 
and containing a term breaking the Uyi(l) symmetry, which reflects the axial anomaly 
in QCD (described by coupling constant qd). When put in a form suitable for the 
bosonization procedure (with four quark interaction only) it results in the following 

§ Actually, in 23] the systematic investigation of the phase diagram of strongly interacting matter as 
a function of temperature and chemical potential was undertaken, which we shall not repeat here in 
the nonextensive version as this would take us outside the limited scope of the present paper. We plan 
to consider it elsewhere. 
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effective lagrangian: 

'Ceff = qin^d, - m) q+Sab [{qX^q) {qX'q)] +Pab [{q^l5X''q) {q^l5X\)] (1) 
where m = diag (m„, nid, m^) and Sab and Pab are projectors, 

Sab = 9sSab + goDabc {qX^q) , (2) 
Pab = gsSab - QoDabc {qX^q) (3) 

with Dabc being the SU(3) structure constants dabc for a,b,c = (1,2,..., 8) whereas 
Doab = —Sab/y/Q and Dqoo = a/2/3. We work with q = {u, d, s) quark fields with three 
flavors, Nf = 3, and three colors, Nc = 3, A" are the Gell-Mann matrices, a = 0, 1, . . . , 8 



and A° = y |l- Integrating over the quark fields in the functional integral with Ceff 
one gets an effective action expressed by the natural degrees of freedom of low energy 
QCD in the mesonic sector, namely a and ip (the notation Tr stands for taking trace 
over indices Nf and Nc and integrating over momentum) : 



Weff[<^,(r] 



(4) 



— zTr In 



Z7^9^ - m + OaX" + {i-i^){^aX'') 



The first variation of Weff leads to the gap equations for the constituent quark 
masses Mi. 

Mi = mi- 2g^ {qiqi) - 2g^ (q^qj) {q^qu) , (5) 
with cyclic permutation of i,j,k = u,d,s and with the quark condensates given by 



(qiqi) = — ^Tr[5'j(p)] {Si{p) is the quark Green function); rrii denotes the current mass of 
quark of flavor i (notice that nonzero go introduces mixing between different flavors). 

Let us consider a system of volume V, temperature T and the i^^ quark chemical 
potential /ij characterized by the baryonic thermodynamic potential of the grand 
canonical ensemble (with quark density equal to pi = Ni/V, the baryonic chemical 
potential ps = ^{fJ'u + ^8 + Ps) and the baryonic matter density as ps = l{pu + Pd + Ps)), 

fi(T, V,p,)=E-TS- J2 /^^^- (6) 



The internal energy, E, the entropy, 5", and the particle number, Ni, are given by 

+ rriiMi 



(here Ei = ^Mf+p^): 



E 



Nr 



p dp- 



Ei 



-(1 - Hi- Hi) 



- 9sV ^ {{qiqi)f -2gDV{uu){dd){ss), 

i=u,d,s 
i=u,d,s 

where 5* = [rii In nj + (1 — Ui) ln(l — nj)]+[raj ^ 1 — 

Nc f 
Ni = -^V / p^dp {rii - Hi) . 



rii 



(7) 
(8) 

(9) 
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The quark and antiquark occupation numbers, rii and rzj, are 

= exp[/3(i?, -/.,)] + !' = exp [{p{E, + f,,)] + 1 

With these occupation numbers one can now calculate values of the quark condensates 
present in Eq. 



1 — Ti,- — n,) 



dp. (11) 



i=u,d,! 

Eqs. ([5]) and (ITT!) form a self consistent set of equations from which one gets the effective 
quark masses Mj and values of the corresponding quark condensates (once a temperature 
and chemical potential are given). 

The values of the pressure, P, and the energy density, e, 

P(,..T) = -^. ,(^„T) = ^ (12) 

are defined such that P(0, 0) = e(0, 0) = 0. 

In order to illustrate the g-dependence of the chiral phase transition for zero 
chemical potential, in the present work we are only concerned with 7r° and a mesons. 
Their effective masses can be obtained from the effective action (j4j) by expanding it over 
meson fields and calculating the respective propagators. In the case of a a meson one 
must also account for its matrix structure in isospin space cf., [23]. And thus the mass 
of the 7r° meson is be determined by the condition: 

Dj{M^o,O) = 0. (13) 

where D~o is the inverse of the meson vr^ propagator, 

DjiP) = l-P^oj!^^iP), (14) 

PttO = gs + 9d {qsqs) (15) 

and operator J^u is given by well defined integrals [27] (cf., also appendix of ref.([23])). 
The procedure for obtaining the effective mass of the a meson is analogous with the 
only difference that solving the condition D~^{M^, 0) = we use approximate form of 
a propagator (see [27] for details). Finally, the model is fixed by the coupling constants, 
gs and go, the current quark masses, m^, and the cutoff A, which is used to regularize 
the momentum space integrals [Jj]. 



3. Nonextensive NJL model - g-NJL 
3.1. Motivation 

As a motivation for study of nonextensive version of the NJL model, the g-NJL 
model, let us notice the following. The NJL model [23] is formulated in the grand 

II For numerical calculations we use the same parameter set as that in |23j : uiu = rud = 5-5 MeV, 
rus 140.7 MeV, gsA^ = 3.67, gcA^ = -12.36 and A = 602.3 MeV. It has been determined by fixing 
the values AU = 135.0 MeV, Mk = 497.7 MeV, U = 92.4 MeV, and M^/ = 960.8 MeV. For the quark 
condensates at T = we obtain: (g„g„> = {qdqd) = -(241.9 MeV)^ and {qsqs) = -(257.7 MeV)^, and 
for the constituent quark masses = Md — 367.7 MeV and Ms = 549.5 MeV. 
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canonical ensemble and assumes the additivity of some thermodynamical properties, 
especially entropy. This is a very strong approximation for the system under the 
phase transition where long range correlations or fluctuations are very important. One 
could, alternatively, consider equilibrium statistics using microcanonical ensembles of 
Hamiltonian systems, whereas canonical ensembles fail in the most interesting, mostly 
inhomogeneous, situations like phase separations or away from the thermodynamic limit 
[3T] . The alternative way to describe the non additivity of interacting systems which 
have long range correlations (including long range microscopic memory) or long range 
microscopic interactions is to use g-statistics [32] . 

Let us illustrate this with two examples. At first notice that in the NJL model 
(which in the mean field approximation is given by the effective lagrangian ([1])) one 
introduces a strong attractive interaction between a quark and antiquark represented 
by couplings gs and go] usually assumed to be independent of the temperature T. 
However, this interaction induces the instability of the Fock vacuum of the massless 
quarks which, in turn, results in the non-perturbative ground state with nonzero {qq) 
condensates and in the breaking of chiral symmetry endowing constituent quarks with 
finite masses. This effect takes place in some range of temperatures so as to control it one 
allows in some cases for a temperature dependent coupling constant gn as, for example, 
|27] . It was assumed there that gn, corresponding to breaking of axial symmetry ?7a(1); 
is given by 



where Tq is a parameter (To = 100 MeV in [27]). As shown in [27], depending on the 
assumed value of Tq, chiral symmetry starts earlier. The retardation effect introduced by 
Eq. (fT6|) violates the simple extensivity of the system, therefore it calls for an effective 
nonextensive description provided by g-statistics [2j. 

The second example concerns the description of the spinodal region in the NJL 
models in which one observes a coexistence of two phases: (a) - the phase with broken 
chiral symmetry and with massive quarks (m ~ 300 MeV) and large negative qq 
condensates which constitute the physical vacuum, it develops for small density; {b) 
- for high density the qq condensates disappear and quarks are almost massless (m ~ 5 
MeV). The highest point on the temperature scale of the coexistence curve, Tcru, is 
the critical point. Of special importance is the fact that, within g-statistic, one can 
discuss the occurrence of negative specific heat in a nonextensive system which has 
an equilibrium second order phase transition [33]. According to this analysis, the 
specific heat is negative in a transient regime and corresponds to meta-stable states. 
Exactly such metastable states are observed in the NJL model during the spinodal 
phase transition below the critical temperature [23] . 

Finally, let us notice that the NJL model does not contain color and therefore 
does not produce confinement. Therefore, resigning from the assumption of additivity 
in this case and introducing a description based on the nonextensive approach, which. 




(16) 
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according to [[8], can be understood as containing some residual interactions between 
considered objects (here quarks) seems to be an interesting and promising possibility. 

3.2. Formulation of the q-NJL 

The nonextensive statistical mechanics proposed by Tsallis [2] generalizes the usual 
BG statistical mechanics in that entropy function (we use convention that Boltzmann 
constant is set equal to unity), 

w w 

SbG = -^PMP^^ Sq = -^pllllqPi, (17) 
1=1 i=l 

Sq Sq=i = Sbg for g 1. Here, q is the nonextensive parameter and In^p = 
jpi-g _ i'^ _ qy "Yhe additivity for two independent subsystems A and B (i.e., such 
that = ■ p^) is now lost and takes the form: 

Sf^^ = S^ + S^ + {l-q)Sp^, (18) 

they are called nonextensive 0. 

The relevant point for further consideration is the g-form of quantum distributions 
for fermions (+1) and bosons (—1), which, following (and [22]) we shall take as: 

^ (19) 



eg(/3(E,-/i,))±l' 
where, for g > 1 considered there, 

{[l + (g-l)x]^ if x>0 
(20) 
[l + (l-g)x]^ if x<0 

and X = P{E — fi). The g < 1 case was not considered in [22] whereas [18] advocated 
use of the usual Tsallis cut prescription in this case, i.e., to allow for a given g < 1 only 
for such values of {E, /i, (3) for which [1 + (1 — q)x] > 0. However, in this case we found 
it more suitable to adopt in this case mirror reflection of Eq. (1201) . i.e., that for g < 1 
one has: 

{[l + (g-l)x]^ if x<0 
(21) 
[l + (l-g)x]i^ if x>0 

This is because only then can one treat consistently on the same footing (and 
for all values of x) quarks and antiquarks, which should show the particle-hole 
symmetry observed in the g-Fermi distribution in plasma containing both particles and 
antiparticles, namely that 

nq{E, /5, /i, g) = 1 - n2.g{-E, P, -/i). (22) 

% It is worth knowing that for subsystems with some special probabihty correlations, it is the BG 
entropy for which extensivity is not valid and is restored only for q ^ 1 (one refers to such systems as 
nonextensive [34j ) . 
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This means, therefore, that in a system containing both particles and antiparticles (as 
in our case) both q and 2 — g occur (or, when expressed by a single q only, that one 
can encounter both q > 1 and g < 1 at the same time). These dual possibilities warn 
us that not only the g > 1 but also q < 1 (or (2 — g) > 1 have physical meaning in 
the systems we are considering. This differs our g-NJL model from the g-version of the 
QHD-I model of [2^. 

Notice that for g ^ 1 one recovers the standard FD distribution, n(/i, T). Actually, 
it is important to realize that for T — * one always gets ng(/i, T) — > n(/i,T), 
irrespectively of the value of g [22j. This means that we can expect any nonextensive 
signature only for high enough temperatures. 

In formulating the g-NJL, in what concerns calculations of the modified FD 
distributions of quarks and antiquaries, we follow essentially the steps undertaken in 
[22] where the g-version of the Walecka model of nuclear matter has been formulated 
and investigated using the cutoff prescription proposed in [18]. This allows a comparison 
of the g-version of both approaches. Our g-NJL model is then obtained by replacing 
the formulas of Section [2] with their g-counterparts in what concerns the form of the FD 
distributions. Additionally, when calculating energies and condensates we follow [3, E] 
and use the g-versions of energies and quark condensates replacing Eqs. ([7]) and (fTT!) 
by: 

i=u,d,s 

- 9sVj2 {^m^)^^ - 2gnV{uu)^{dd)^{ss)^, (23) 

i=u,d,s 

and 

i=u,d,s 

On the other hand, again following [9l[5], densities which are given by the the g-version 
of Eq. ([9]) are calculated with n^'s (not with n^, as in (l23l) and in ( 12^ ). The pressure 
for given g is calculated using the above Eg and the g-entropy version of Eq. ([8]) with 

(cf. m) 

Sg=[nlilnqngi + {I - n,i)'^ln,(l - Ugi)] + {ugi l-Ugi} . (25) 

4. Results 

In our study we concentrate on two features of the g-NJL model, namely chiral 
symmetry restoration and spinodal decomposition, the results for which we discus in 
what follows. Because our goal was to demonstrate the sensitivity to the nonextensive 
effects represented by the |g — 1| 7^ 0, we do not reproduce here the whole wealth of 
results provided in [23], but concentrate on the most representative. 



p dp 



p^ + rriiMi 



n 



p'M, 



n„ 



qiJ 



dp. 



(24) 
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4.1. Chiral symmetry restoration in the q-NJL 

Chiral symmetry restoration is best illustrated by the temperature dependence of the 
quark condensates and effective masses of quarks presented in Fig. [T]for different values 
of parameter q. In addition, in Fig. [2] the masses of vr and a mesons for different q are 
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Figure 1. (a) - quark condensates and (b) - effective quark; masses as functions of tlie 
temperature for different values of the nonextensive parameter q (q=l correspond to 
Boltzmann Gibbs statistics). 

presented as function of temperature tEI- They were calculated assuming zero chemical 
potentials and solving numerically the g- version of gap equations and flTT]) . There 
is a noticeable difference for q < 1 and q > 1 cases: whereas g < 1 leads to chiral 
symmetry restoration starting earlier but in general following the usual shape, for g > 1 
it is smeared, starting earlier and ending later. The effects caused by nonextensivity are 
practically invisible for heavier quarks. Also, as seen in Fig. 2, the degeneracies of vr 
and a mesons starts earlier for g < 1 and much later for g > 1. Here we derive the mass 
spectra with the help of the g- version of Eq. (fT3l) solved numerically. The temperature 
for which the a mass reaches a minimum is shifted to smaller values for g < 1 and to 
larger ones for g > 1 by amount depending on the value of g — 1. Also the final value 
of the masses is larger for g < 1 and smaller for g > 1, the actual amount depends on 
the value of |g — 1|. When interpreting g > 1 as a measure of temperature fluctuations 
[H [m [15] , this would mean that fluctuations dilute the region where the chiral phase 
transition takes place; it is especially visible for the a meson which is still not saturated 
at the temperature T = 50 MeV, see Fig. [2l On the other hand, correlations, which 
according to [161 IE] result in g < 1, only shift the condensates, quark masses and 
meson masses towards smaller temperatures. In our case the supposed fluctuations and 
correlations refer to quarks, not hadrons. It is interesting therefore to note that in [22], 

There is still an ongoing discussion on the meaning of the temperature in nonextensive systems. 
However, in our case the small values of the parameter q deduced from data allow us to argue that, 
to first approximation, Tq ^ T used here and in [22j . In high energy physics it is just the hadronizing 
temperature (and instead of the state of equilibrium one deals there with some kind of stationary state). 
For a thorough discussion of the temperature of nonextensive systems, see |35j . 
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Figure 2. Masses of tt and a mesons as functions of tlie temperature for different 
values of the nonextensive parameter q (q=l correspond to BG statistics). 

where degrees of freedom used are nucleons, a similar shift towards smaller temperatures 
occurs for g > 1, i.e., for temperature fluctuations of nucleons (g < 1 is not considered 
there). Notice that, as mentioned in Sec. [H we are not specifying here what are the 
actual dynamical mechanisms behind such fluctuations/ correlations, we just model them 
by the parameter q. An example of such dynamical effect, the temperature dependence 
of the the respective coupling constants, is mentioned above, cf. Eq. (|T6l) . The same 
remark also applies to the discussion that follows. 

4-2. Influence of nonextensivity on the spinodal decomposition in the q-NJL 

The next point we shall address is the influence of g-statistics on the spinodal phase 
transition discussed in [23]. To this end we must proceed to finite density calculations. 
As in |23], we assume chemical equilibrium in the form of fJ^u = fJ'd = fJ's = /W- This 
allows us to work with nonzero Ni in Eq. ([9]) and thus with nonzero baryon density 
defined as p = | J2i ^i/'^ ■ 

The spinodal phase transition occurs, in general, for finite densities and for 
temperatures below a critical temperature T^t [Ml [I]- Above it we do not observe 
phase transition of the first order but rather a smooth crossover. Below it, for some 
range of densities, we have a region of mixed phases of hadronic and quark matter [23] 
(understood here as phases with very small (current) and very large (constituent) quark 
masses). The first observation is that details of the spinodal phase transition are very 
sensitive to g — 1, much more than it was observed in the previous case. It is best seen 
in the g-dependence of Tcr shown in Fig. [3l it changes by ~ 10 MeV between g = 0.98 
and g = 1.02 used here. The general observation is that for g < 1 the pressure decreases 
and energy increases, whereas for g > 1 one observes the opposite tendency, cf. Figs. H] 
- [3 Such behavior is a direct consequence of the nonextensivity and arises from the the 
(g — 1) term in the nonextensive entropy functional f|T8|) . 

In Fig. Hlwe show the pressure at the critical temperatures Tcr for different values 
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Figure 3. Critical temperature Tcr as function of the nonextensivity parameter q (in 
the range of q considered here). 




p/po 

Figure 4. The pressure at critical temperature Tcr as a function of compression p/po 
calculated for different values of the nonextensivity parameter q. The dots indicate 
positions of the inflection points for which first derivative of pressure in compression 
vanishes. As in [23j for q = 1 the corresponding compression is p/po = 1-67 (and 
this leads to p = 318 MeV); it remains the same for g > 1 considered here (but now 
p = 321 MeV for q = 1.01 and p = 325 MeV for q = 1.02) whereas it is shifted to 
p/p = 1.72 for q < I (resulting in p = 313 MeV for q = 0.99 and p = 307 MeV for 
q = 0.98). 

of g as a function of the compression p/po (with po = 0.17 fm^^). Notice that the 
effect is stronger for g < 1 and that, essentially, the saddle point remains at the same 
value of compression. When one moves away from the critical temperature, one gets 
the pressure versus compression as presented in Fig. [5] for two different temperatures, 
T = 30 and 50 MeV. Notice the occurrence of the typical spinodal structure, which is 
more pronounced for lower temperatures, T = 30 MeV, whereas its sensitivity to the q 
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Figure 5. The pressure calculated for different values of the nonextensivity parameter 
q for temperatures T — 30 Mev (a) and T = 50 MeV (b) as function of the compression 
p/ Po- The curves for q for which the temperature considered is the critical temperature 
are also shown, they correspond to q = 1.19 for T = 30 MeV and q = 1.063 for T = 50 
MeV. 
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Figure 6. Pressure corresponding to the local minimum in Figs. [4] and [5] versus 
temperature T calculated for some selected values of the nonextensivity parameter q. 
The curves end at the critical points. 

parameter gets stronger with increasing temperature. Shown are also curves for g's for 
which the temperature considered coincides with the critical temperature for this value 
of q. The are, respectively, q = 1.19 for T = 30 MeV and q = 1.063 for T = 50 MeV, 
i.e., the corresponding values of q decrease with temperature, as expected from Fig. [3l 
It means that for each temperature (even for very small one) a g > 1 exists for which 
there is no more mixed phase and for which spinodal effect vanishes. This seems to be 
a quite natural effect in the scenario in which g > 1 is attributed to the fluctuations of 
the temperature in a system considered as proposed in [HI H]. On the contrary, effects 
like correlations or limitations of the phase space considered in [T6l [T7] work towards an 
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Figure 7. The energy per particle at the temperatures T = 30 and T = 50 MeV 
as a function of compression p/po calculated for different values of the nonextensivity 
parameter q. 



increase of the Tcr and make therefore the spinodal effect more pronounced. 

Fig. [6] shows the temperature behavior of the pressure Pg defined as the pressure 
at the local minimum in the spinodal curve. The characteristic features are that the 
temperature at which Pg starts to be positive and is shifted towards smaller values with 
increasing q and that it gets deeper into negative values with decreasing q. The effect of 
negative pressure can be best understood invoking the bag model picture of the nucleon 
[29] . Generally, in the phase of hadron gas we observe a decrease of the pressure and the 
critical temperature Tcr with the increase of fluctuations given by q. That phenomenon 
resembles to some extent the behavior of the bag constant for nucleon in the medium 
where the bag pressure decreases with the increase of the chemical potential /i in order 
to get a proper equation of state [37] . In the nuclear thermodynamical models this bag 
constant is modified because the vacuum, in which hadrons are embedded, is modified 
by the residual interaction present in the nuclear medium (acting towards the Wigner 
realization of the chiral symmetry in which masses of tt and a are degenerated). Here 
such a density dependence corrections are introduced by nonextensive effects inside the 
nuclear medium. In that way the nuclear vacuum for the temperatures below the critical 
temperature and critical densities, the usual area of the spinodal phase transition, can 
be properly described effectively by the nonextensive statistics. 

Fig. [7] shows the energy per particle, E/A (cf. Eq. fl23|) for different temperatures 
and for different values of the nonextensivity parameter q. Notice that whereas for 
g < 1 this energy exceeds the usual one (i.e., for g = 1) it gets smaller for g > 1 (this 
effect is especially seen for compressions smaller than two). This is an opposite trend 
to that observed for the corresponding behavior of the pressure. Finally, it is worth to 
observe that the absolute minimum of energy for given temperature does not depend 
on the nonextensivity parameter q and for T = 30 MeV it is located at p/po = 2.45. 
It turns out that to obtain stable state here, i.e., P = 0, one has to choose q = 0.97. 
In such a way the final droplets of quarks |38] in the mixed phase can appear at finite 
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Figure 8. Phase diagram in the g-NJL model in T — /i plane for values of q considered 
before: q — 0.98, 1.0 , 1.02. Solid and dashed lines denote, respectively, the first order 
and crossover phase transitions. The results are presented for three different values 
of the nonextensivity parameter q with the vicinity of the (g-dependent) critical end 
points (CEP) enlarged in the inlet. The crossover phase transition for q — 0.98 and 
for /i ^ takes place for smaller temperature T. 



temperatures. 

Let us close with few remarks. Nonextensive dynamics enter the NJL calculations 
through the quark (antiquark) number distribution functions n^j {figi). These functions 
are connected with the respective quark (antiquarks) spectral functions in NJL model. 
However, deviations from the exponential shape of g-exponents, as defined in Eqs. ( l20l) 
and (12T|) . are negligible for values of q close to unity (in our case 0.98 < q < 1.02). 
It is also important to notice that Eqs. (!20!) and (12T!) are symmetric for q and 1 — q. 
The differences between g < 1 and q > 1 cases observed in our results are then due 
to our way of defining the energy fl23|) and entropy fl25l) . which, following [9l E], we do 
by using n'^gj now for g < 1 distributions 

n'^gi and n'^gi are closer to unity than Ugi and figi, therefore the absolute values of quark 
condensates (as given by Eq. ( !24l) ) begin to decrease for q = 0.98 at lower temperature 
in comparison with the g = 1 case, see Fig. [H The corresponding energy, see Fig. 
[71 is bigger, which means that g < 1 introduces some residual attractive correlations 
which rise the energy and lead to hadronization occurring at lower temperature (see 
Figs. [Hand [2]). On the other hand, g > 1 introduces fluctuations which decrease the 
effective occupations {n'^gi and n'^gi) and the energy, and smears out the chiral phase 
transition, see Fig. [2j In Fig. [8] we present our phase diagram in the /i — T plane for 
different nonextensivity parameters considered here with positions of the corresponding 
critical end points (CEP) for different values of g clearly indicated. The overlap of 

* It is worth to notice that in [22], which considers only the q > 1 case and uses number distributions 
without powers of the significant effects were obtained only for much bigger values of the nonextensive 
parameter q = 1.2. 
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curves observed in Fig. [H] (inlet) indicates how critical end point is smeared to a kind 
of critical area. This is because fireballs created in different events can have different 
values of q (representing, as mentioned before, action of all factors responsible for the 
departure of our system from the usual BG approach - not specified here in detail but, 
in general, resulting in specific correlations of quarks or fluctuations of temperature 
mentioned before). Therefore when analyzing experimental data one most probably 
will encounter such a critical area instead of well defined CEP. 

5. Summary and conclusions 

To summarize: we have investigated the sensitivity of the mean field theory of the NJL 
type presented in |23j to the departure from the conditions required by the application of 
the BG approach. To this end we have used the Tsallis version of nonextensive statistical 
mechanics [2] with a new parameter q, the phenomenological nonextensivity parameter, 
such that q—1 quantifies departure from the BG situation (which is recovered for g ^ 1). 
As result we have obtained a nonextensive, g-NJL model with q being a new parameter 
summarizing action of a number of yet undisclosed factors, which should be fitted to 
the data. Our investigation was prompted by recent investigations of similar effects in 
quantum hadrodynamics [22]. In fact, we have used the same nonextensive version of 
the standard FD distributions (discussed in [H]), but this time applied to quarks and 
antiquarks and for both g > 1 and g < 1 cases. On the other hand, when calculating 
energies, quark condensates and densities we followed prescription advocated in |9l |5] . 

We have investigated two possible scenarios corresponding to g > 1 and g < 1, 
respectively, which, as mentioned, correspond to different physical interpretations of the 
nonextensivity parameter. For g < 1 (usually connected with some specific correlations 
[16] or with fractal character of the phase space [T7]) we observe decreasing of pressure, 
which reaches negative values for a broad (g-dependent) range of temperatures and 
increasing of the critical temperature % The g > 1 case (usually connected with 
some specific nonstatistical fluctuations existing in the system [HI H] [|J we observe 
a decreasing of the critical temperature, Tcru (cf. Fig. [3]), and therefore in the limit of 
large g we do not have a mixed phase but rather a quark gas in the deconfined phase 
above the critical line (on the contrary, the compression at critical temperature does 
not depend on g (cf. Fig. HI). As in [22] the resulting equation of state is stiffer (in the 
sense that for a given density we get bigger pressure with increasing g). As expected, 
the observed effects depend on the temperature and tend to vanish when temperature 
approaches zero. Finally, as shown in Fig. [8|, the nonequilibrium statistics dilutes the 
border between the crossover and the first order transition. 

We would hke to end by stressing that our results could be of interest for 
investigations aimed at finding the critical point in high energy heavy ion collisions 

tt It acts therefore in the same way as including of the Polyakov foop into the NJL model [39| . 

ft Actually, in 7J it was shown that if one could find a dynamical source of such fluctuations and 

introduce it to the model then the new q would diminish considerably, eventually becoming unity. 



Nonextensive effects in the Namhu-Jona-Lasinio model of QCD 



16 



[Ij or when studying particularities of the equation of state (EoS) of compact stars [ID] . 
The fact that they do depend on the parameter q used means that the exact position of 
such a point or the type of the or the shape of EoS could be quite different from that 
naively expected. 
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